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Abstract. The observations of high energy 7-ray emission from the Galactic center 
(GC) by HESS, and recently by Fermi, suggest the cosmic ray acceleration in the GC 
and possibly around the supermassive black hole. In this work we propose a lepton- 
hadron hybrid model to explain simultaneously the GeV-TeV 7-ray emission. Both 
electrons and hadronic cosmic rays were accelerated during the past activity of the 
GC. Then these particles would diffuse outwards and interact with the interstellar gas 
and background radiation field. The collisions between hadronic cosmic rays with gas 
is responsible to the TeV 7-ray emission detected by HESS. With fast cooling in the 
strong radiation field, the electrons would cool down and radiate GeV photons through 
inverse Compton scattering off the soft background photons. This scenario provides a 
natural explanation of the observed GeV- TeV spectral shape of 7-rays. 
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1. Introduction 

It is well known that the Galactic center (GC) region has a very complex astrophysical 
enviroment and is rich in various kinds of objects, and is a good library for the study of 
astrophysical phonomena. A supermassive black hole with mass ~ 4 x 10 6 M Q , Sgr A*, 
lies in the GC. Although the GC is rather quiet nowadays, frequent flares were observed 
in the X-ray as well as near infrared (NIR) bands [H [2j [3], which means the existence 
of continuous weak activities of the black hole. 

As a consequence high energy particles could be accelerated during such kinds 
of activities and imprinted in the 7-ray sky. The very high energy (VHE) 7-rays 
from the GC have been observed by several atmospheric Cerenkov telescopes such as 
CANGAROO [4j, VERITAS [5], HESS [8], and MAGIC [9J. Recently GeV 7- 
ray emission from the GC was also revealed by the spatial detector Fermi Large Area 
Telescope (Fermi-LAT, [TO])- The observations in both TeV and GeV bands showed 
that the 7-ray emission tended to be stable without significant variability [El El [10] • 

Serveral models have been proposed to explain the 7-ray emission observed at the 
GC, including hadronic models [TU [El CEJ [101 HU HH] and leptonic ones [HI [T7J. An 
issue prevents a simple explanation to the GeV- TeV spectra of the GC is that the GeV 
7-ray emission is not consistent with the direct extrapolation of TeV 7-rays to the low 
energy range. Chernyakova et al. [TO] proposed that the strong energy dependence of 
the diffusion coefficient would seperate the low energy and high energy particles into 
two different diffusion regimes and the GeV- TeV spectral shape of the 7-rays could be 
reproduced. In [15] the stochastic acceleration of particles in a two-phase interstellar 
medium (ISM) was employed to generate two distinct populations of protons to explain 
the data. In [17] the Fermi-LAT 7-ray emission was suggested to originate from a 
population of electrons ICS off the soft background photons and the VHE 7-ray emission 
was thought to be from different sources instead of Sgr A*. 

The morphology study of the VHE 7-ray source, HESS J1745-290, by HESS showed 
no spatial extension and the upper limit of the angular size of ~ 1'.3 was set [18] . Note 
there is another diffuse component of the 7-ray emission along the ridge [HI [20] , which 
is not discussed in this work. It may imply that the high energy particles are confined 
in several parsec regions around the GC [TO1 [15] or the high concentration of the target 
gas in the inner region [14] . 

In this work, we propose a hybrid model to explain the GeV and TeV 7-ray emission. 
It is natural to expect that both protons and electrons can be accelerated during the 
GC activity. Therefore without finely tuning the environmental parameters, we should 
expect the existence of two components of the 7-ray emission, which is just the case 
shown by the GeV- TeV observations. A simple expectation is that high energy protons 
interact with the ISM is responsible to the VHE 7-rays, and the electrons, which may 
cool down in the GC radiation field and/or magnetic field, produce the GeV 7-rays 
through inverse Compton scattering (ICS) off the background photons. In Sec. 2, we 
present the detailed modeling of this picture. Sec. 3 is the discussion and the main 
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conclusions are summarized in Sec. 4. 



2. Model and results 



Although the overall behavior of Sgr A* is quite silent, the observations in X-ray 
and infrared bands indicate that it still has continuous weak activities [2lT [221 [U [2]. 
During such kinds of avtivities, the accretion of stars and gas by the supermassive black 
hole could be effictive to accelerate particles. We assume that both the protons and 
electrons were accelerated instantaneously during such flare events. Then these particles 
would diffuse away from the acceleration site, and radiate during the propagation. The 
propagation equation for both protons and electrons can be written as 

where 0(r, E, t) is the propagated flux as a function of space, energy and time, D(E) is 
the diffusion coefficient, N(E) is the injection spectrum of particles at t — and r = 0, 
and b(E) = dE/dt is the energy loss rate, which is important for electrons but negligible 
for protons. The energy dependence of diffusion coefficient is assumed to be power-law 
of rigidity R, D(R(E)) = (3D (R/4: GV) S , in which 6 is the spectral index and j3 is the 
particle velocity. In this work we assume 5 = 0.3, which is close to a Kolmogorov type 
of the ISM turbulence. 

The solution of the above equation for electrons is [23] 

6(rEt)- N <W b W exp (2) 

where Ei is the initial energy of the electron whose energy is cooled down to E at time 
t, and r<iif(E, t) = 2\/ Au, in which 

is the effective diffusion length when the particle energy decreases from Ei to E. 

For protons we neglect the energy loss term, and the solution of the propagation 
equation is even simpler 



N p (E) ( r 



,.2 



E, t) = - ^_ [ 3 exp ( -— ) , (4) 



where a(E,t) = \j2D(E)t is the effective diffusion length within t. 

The energy loss of electrons includes ICS off the soft background photons, 
synchrotron radiation in the magnetic field and bremsstrahlung radiation in the ISM. 
The soft photons consist with infrared from dust and optical from stars. In [22J the 
measurement of infrared emission in the central 1.2 pc was presented. According to the 
data, the expected optical emission was proposed based on assumptions of the optical 
photon energy density [23J [161 EZ]- Taking the uncertainties of the expected optical 
emission into account, we adopt the two models given in [211 E] i n this work. The 
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Figure 1. Infrared and optical photon energy density spectrum in ~ 1 pc of the GC 
region. The solid line is adopted from [33] which is consistent with the observational 
data [22] . The dashed and dash-dotted lines are the model expected results of optical 
emission presented in [23] and [T7] respectively. This figure is a reproduction of Fig. 1 
in [17]. 



intensity of the soft photons is shown in Figure [T] The magnetic field strength varies 
with the radius away from the GC [25]. For the very small region (r ~ pc) the average 
magnetic field strength is adopted to be ~ 10 2 /iG [25] H7) . As for the ISM density, 
there is large uncertainty in the GC region. Typically the adopted ISM density within 
~pc region is of the order 10 3 cm -3 (Till E]. Compared with the ICS energy loss in the 
strong radiation field, the bremsstrahlung energy loss is negligible |26j . 

We calculate the cooling time of the electrons, r = E/(dE/dt), in the above 
background photons field and magnetic field, as shown in Figure |2] Here the optical 
emission is adopted to be the one of Kusunose & Takahara (2012) [IT]. It is shown that 
the energy loss of electrons is dominated by the ICS up to > 10 TeV. 

In Figure [3] we show the integral spectra of both electrons and protons after time 
t of the injection. The integral radius is taken to be 1.2 pc for electrons and 3 pc for 
protons respectively. For electrons we integrate the particles within radius ~ 1.2 pc since 
the high background photon density is expected to exist in that region, and most of the 
radiation should come from such a region. However, for protons, we adopt a some larger 
integral radius, ~ 3 pc, according to the spatial extention upper limit of HESS J1745-290 
derived by HESS [IS] . The diffusion coefficient is Dq = 10 27 cm 2 s _1 , and the injection 
spectra of both electrons and protons are adopted as E~ 19 . We can see from this figure 
that the cooling of electrons will result in a cutoff of the spectrum, above which the 
electrons can not survive from the energy loss. The cutoff energy depends on the time. 
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Figure 2. The cooling time of electrons in the background photon field model of 
Kusunose & Takahara (2012) [17 and magnetic field with B = lOO^G. The ISM 
density is adopted to be 10 3 cm~ 3 to calculate the bremsstrahlung energy loss. 



10° br 



c 
z> 

£• 

CO 

I— 

-t— ; 

< 

LU 



10" 



LU 



10" 




100 yrs(Lepton) 
200 yrs(Lepton) 
300 yrs(Lepton) 
100 yrs(Proton) 
200 yrs(Proton) 
300 yrs(Proton) 



10 L 



10' 



10' 



10" 



E (GeV) 



Figure 3. Propagated electron (solid) and proton (dashed) spectra for observational 
time t = 100, 200 and 300 (from top to bottom) years respectively. The energy loss of 
electrons adopted in this calculation is shown in Figure [2] 



The younger the injection is, the higher energy of electrons we can observe. The spectral 
shape also has an evolution with the time. For shorter time, the relative fraction of low 
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energy particles is smaller due to slower propagation of low energy particles, which 
results in harder spectra. For protons there is also a time evolution of the spectra. It 
reflects the fact that for longer time, more high energy particles will diffuse out of the 
integral region and result in a soft spectrum. 

Then we discuss the 7-ray emission of those electrons and protons in the background 
radiation field and the interstellar medium. For electrons we consider only the ICS 
emission and neglect the bremsstrahlung radiation. The Klein-Nishina cross section of 
ICS is employed. For the 7-ray emission from pp collisions, we use the parameterization 
given in [27J. 
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Figure 4. Calculated 7-ray spectra of the hybrid model, compared with the 
observational data by Fermi-LAT [TU] and HESS [28 . Different lines show the effect of 
observational time t (1 00, 200 and 300 yr from top to bottom). The diffusion coefficient 
is Do = 10 27 cm 2 s _1 . The soft background photon model is Kusunose & Takahara 
(20f2) [17]. 

Figure IHpresents the results of the calculated 7-ray spectra and the comparison with 
the observational data [10l|28]. The injection spectra are adopted to be oc E' 1 - 9 , and the 
diffusion coefficient is adopted to be D = 10 27 cm 2 s -1 . The background photon model 
is adopted to be Kusunose & Takahara (2012) [T7j. Three different observational time 
t = 100, 200 and 300 yr are shown. The observational time determines the cutoff energy 
of the electron spectrum (see Figure [3]). We find that for the Kusunose & Takahara 
(2012) photon field t ~ 200 yr can match the data welQ. Another effect of the time t, 
is the overall normalization. If t is larger, more particles diffuse out of the ~pc region, 

I Here we do not intend to fit the data to find the best parameters, but just to show the rough 
comparison between the model and data. Except the first data point of Fermi-LAT, we actually find 
rather good description to the data of the model (dashed line) . There could be systematic uncertainty of 
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and we expect a decrease of the overall normalization of the integral particle spectra. 
For t = 200 yr we find the total energy of electrons above 1 GeV is ~ 5 x 10 47 erg to 
match the Fermi-LAT data. The electron-to-proton ratio K ep is about 0.22(n/10 3 cm -3 ) 
in this case. Note, however, a significant part of high energy protons (E > 10 TeV) will 
diffuse out of the 3 pc region in this case, as can be seen from Figure El Therefore to be 
consistent with the morphology of the source, we may require the gas density is higher 
in inner 3 pc region than outside, as the scenario of [14J. 




E (GeV) 

Figure 5. Same as Figure|4]but for smaller diffusion coefficient Do = f0 26 cm 2 s _1 . 

To be consistent with the cutoff behavior of VHE 7-rays observed by HESS, we 
employ an exponential cutoff of the proton spectrum with cutoff energy E c = 200 TeV. 
It is possible that the pair production of VHE photons in the background radiation 
field could also be responsible to the cutoff of the 7-ray spectra [29], [30] • However, our 
calculation shows that for the adopted soft background photon field in ~pc region, the 
pair production attenuation is negligible. Thus the cutoff should be understood as the 
acceleration limit of the flaring event. The maximum energy protons can achieve for 
diffusive shock acceleration is [31] 

£ max ~ eBR « 10 21 (2) (^) eV > (5) 

where B is the magnetic field and R is the size of the acceleration region. As in [3T] , 
we assume the acceleration takes place within 10 Schwarzschild radii (R g ~ 10 12 cm) of 
the black hole. To accelerate protons to above 100 TeV requires a magnetic field ~G 



the data analysis of Fermi-LAT, especially at low energies, and adjusting properly the model parameters 
may also further improve the fit. 
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in the acceleration region. Such a condition could be reached in the very central region 
of the GC [31] . On the other hand, if the acceleration takes place in larger regions, the 
required magnetic field could be smaller. 

Figure [5] shows the results for even smaller diffusion coefficient D Q = 10 26 cm 2 s -1 . 
In this case most of the protons up to 100 TeV may still be confined in several pc 
region around the GC. As a simple estimate, the diffusion length for protons is about 



similar with the injection spectrum. Therefore we adopt the injection spectrum oc E~ 2,2 
to be consistent with the HESS data. The results show a good fit to the GeV-TeV 
observational data for t = 200 yr too. It is shown that for protons different time t 
leads to small differences of the results, which means that indeed most of the protons 
are confined. We also note that the differences of the ICS spectra with respect to t 
are larger than that of 7r° decay spectra. This is simply because the integral radius of 
electrons is smaller than that of protons. The total energy of electrons above 1 GeV is 
~ 8 x 10 46 erg in this case, and the electron-to-proton ratio K ep is 0.05(n/10 3 cm -3 ). 

Finally, the results for Mezger et al. (1996) soft background photon model [21] are 
shown in Figure [6l Other conditions are the same as Figure HI Compared with Kusunose 
& Takahara (2012) photon field, the optical emission is stronger and the cooling time 
for electrons is shorter. Thus for t ~ 150 yr the cutoff energy of electrons is proper 
to fit the Fermi-LAT data. However, the energy spectrum of ICS emission does not 
differ significantly from that shown in Figure HI This is because in stronger background 
radiation field the cooling is more significant, and a larger cooling is just canceled out 
when calculating the ICS 7-ray spectrum in such a radiation field. The total energy of 
electrons above 1 GeV is ~ 1.6 x 10 47 erg, and the electron-to-proton ratio K ep is about 
0.11(n/10 3 cm" 3 ). 

3. Discussion 

The diffusion coefficient adopted in this work is Dq ~ 10 26 — 10 27 cm 2 s -1 , which is 
much smaller than that in the disk and halo (~ 10 28 — 10 29 as induced from the cosmic 
ray transportation). This might be due to the higher magnetic field and more turbulent 
ISM in the GC region. The lower limit of the diffusion coefficient should be the Bohm 
limit, which is D Bohm = pu/3eB ~ 3 x 10 21 (p /TeV) (B / fiG)- 1 cm 2 s" 1 [32j [331 [33] . The 
diffusion coefficient adopted in this work is well larger than the Bohm limit in the whole 
energy range discussed. 

In the calculation of proton propagation, the energy loss of protons via interaction 
with the ISM is neglected. As a rough estimate, the average collision probability of 
one projectile proton is about navr ~ 0.01(ra/10 3 cm" 3 )(r/200yr). Thus for the typical 
parameters adopted in this work, the collisional energy loss of protons is indeed negilible. 

The secondary electrons/positrons from charged pion decay due to pp collision 
could also contribute to the 7-ray emission. As a very rough estimate, for one pp 
collision, 7-ray photons will take 1/3 of the energy of all the pions, and e will take 




The accumulative spectrum of protons is 
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Figure 6. Same as Figure @] but for the Mezger et al. (1996) soft background photon 
model [21]. 

2/3x1/3 = 2/9 of the energy (assuming energy equipartition between e ± and neutrinos) . 
Therefore the total energy of secondary e should be lower than that of n° decaying 
7-rays. Furthermore, the soft radiation field is restricted in 1.2 pc region around the 
black hole, while the pp collision occurs in ~ 3 pc region. Therefore the contribution of 
the secondary to the 7-rays through ICS emission should be much smaller than that 
of neutral pion decay. 

The electron component will also produce synchrotron radiation in the magnetic 
field, and there might be constraints from the radio to X-ray data. As shown in [17] the 
synchrotron emission in ~ 100 /zG magnetic field may exceed the radio measurement of 
Sgr A* for the quiescent state [35]. However, the spatial scale of the present study is 
arc-minutes (pc), instead of arc-seconds as the radio observations show. The expected 
synchrotron emission is consistent with the data-based result of pc scale radio emission 
as given in [24J. 

Early in 2012, one bright flare of GC is observed by HETGS onboard of the Chandra 
X-ray observatory [3]. The total energy in 2 — 10 keV band was approximately 10 39 erg 
[3]. Assuming the total energy of this accretion event is about 4 orders of magnitude 
higher [24], ~ 10 43 erg, and ~ 10% of it converts to acceleration of cosmic rays, we 
give the expected 7-ray emission of such an event within the current framework of the 
model, as shown in Figure [7J The solid lines are same as that in Figure HJ The dashed 
line is the expected flux of the ICS emission, with the same parameters of the diffusion 
coefficients, injection spectral index and electron-to-proton ratio as that of Figure HI 
The 7-ray emission at this stage is dominated by the ICS emission and the hadronic 
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component is much lower. We can see that the total flux of 7-rays is too low to be able 
to be detected by the current VHE 7-ray detectors. We should keep in mind that this 
estimate is very rough and suffers from large uncertainties of the assumption of total 
energy output of the flare and the energy fraction goes into particle acceleration. This 
gives us an impression that how large a flare is needed in order to produce the observed 
7-rays. 




E (GeV) 

Figure 7. The expected 7-ray emission for the flare event of 2012. We assume the 
total power of the flare event is 4 orders of magnitide high than the 2—10 keV energy. 

4. Conclusion 

In this work, a hybrid model is proposed to explain the GeV 7-ray emission observed by 
Fermi-LAT and TeV 7-ray emission observed by HESS telescope, of the GC. The current 
7-ray emission could be originated from one flare with total energy > 10 48 ergs at ~ 10 2 
years ago. Both the protons and electrons were accelerated to very high energies in this 
flaring event. Furthermore the flaring event produced a thermal bath of soft photons 
in the GC region, with typical scale ~pc. High energy particles could then diffuse out 
of the acceleration site and interact with the ambient medium and radiation fiels. The 
hadronic collisions between protons and ISM gives rise to the TeV 7-rays. At the same 
time, electrons would lose energy through ICS and synchrotron radiation. The ICS 
photons could be responsible for the GeV 7-ray emission observed by Fermi-LAT. 
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